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Abstract 

Using the method of the Laplace transform, we consider fractional oscillations. They 
are obtained by the time-clock randomization of ordinary harmonic vibrations. In 
contrast to sine and cosine, the functions describing the fractional oscillations exhibit 
a finite number of damped oscillations with an algebraic decay. Their fractional 
differential equation is derived. 
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1 Introduction 



According to fractional calculus tlie liarmonic oscillator is a particular case of 
the fractional oscillator [Tp|3] . It suffices to say that the fractional oscillator 
equation is a generalization of the classical harmonic equation by replacing 
the second-order derivative by a fractional one. If the harmonic oscillator is 
described by elementary harmonic functions, the fractional oscillator solution 
is expressed in terms of Mittag-Leffler functions. The interrelation between 
these oscillators has clarified with the consideration of temporal subordina- 
tion. Then the development of an oscillatory process is governed by its own 
internal clock (operational time) that is not synchronized with physical (de- 
terministic) clock [4J. As has been shown in [5J that the fractional oscillator 
results from an ensemble average of harmonic oscillations under stochastic 
time arrow. In fact, the fractional oscillator can be considered as a time-clock 
randomization (subordination) of the conventional harmonic oscillator. The 
new time clock is the continuous limit of the discrete counting process, when 
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the time variable is a sum of random temporal intervals belonging to an a- 
stable distribution. Hence there exists a direct correspondence between the 
functions of both oscillators. 



In particular, the cosine function cos(t) is connected with the one-parameter 
Mittag-Lefiler function Ea{—t"')- Its features has been studied in detail (see, 
for example, [1] and references therein). Passing to the limit a — > 2, this 
Mittag-LefHer function transforms to cos(t). One of very interesting proper- 
ties of the one-parameter Mittag-LefHer function is a finite number of zeros. 
Therefore, the function has a finite number of damped oscillations with an 
algebraic decay. For completeness it is necessary to investigate also the "sine" 
correspondence of fractional oscillations. The function is briefly mentioned 
in [3j without any detailed analysis. In this connection it should be pointed 
out that the displacement of harmonic oscillator and its momentum are writ- 
ten by means of ordinary harmonic functions. Similar values (displacement and 
momentum) may be formulated for fractional oscillator too. Their consider- 
ation is of great interest. They are just expressed in terms of Mittag-Leffler 
functions mentioned above. The aim of this paper is to study their features in 
the context of fractional oscillator. 

The paper is organized as follows. In Section |2] we start our analysis with the 
time-clock randomization of the ordinary harmonic oscillator. This allows us to 
use the integral relation between the harmonic functions and the Mittag-Leffler 
functions of the fractional oscillator. Our main interest will be focused on the 
Mittag-Leffler function having the sine function as a limit. Next we investigate 
its properties in Section [3] and derive a fractional differential equation for 
which this function is its solution (Section H]). As a comparison, the "cosine" 
correspondence of fractional oscillations will be given simultaneously. Section [5] 
is devoted to a consideration of fractional oscillation zeros. 



2 Probabilistic point of view 

If one randomizes the time clock in accordance with [S] , the fractional oscilla- 
tions are written by means of 
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/(t,r) coscurdr = E^^i{-ujH'^) = E„(-a;V) , 
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where < a < 2, 

oo 

is the Mittag-Leffler function, T{z) the gamma function, and 

Br 

determines the probabihty to be at the internal time r on the real time t. 
Here Br denotes the Bromwich path (a line Res = a with a value cr > 1), 
and j = 

Since the function A{t) has been reviewed in [1], we devote our attention more 
to the function B{t). The latter has a Laplace inversion integral 



B{t) = — / e^*^ -ds. 

Br 



The fractional oscillations exist for 1 < a < 2. It is that the case is of interest 
to us and will be discussed below. 

For convenience one puts co = 1, then one denotes B(t) as ia(t). It should 
be pointed out that for w = 1 the function A{t) is equal to ea(t) accord- 
ing to the notations of [Ij. As the Laplace image of the function ea{t) = 
X]fcLo(~^°)^/r(tt^ + 1) is + s"), it is easy to determine the relation 

between ea{t) and ia{t) via the fractional integral of a/2 order, namely 

t 



= -r^'e^it) = —— / (t - r)-i+"/2 e„(r) dr . 
T{a/2) 







Recall that for a = 1,2 the function ea{t) takes the form 
ei(t) = e~*, e2(t) = cost . 



Define the same in reference to the function ia{t) for a = 1,2. At once it 
is obvious that i2{t) = sint. Using the Efros' theorem, we find ii{t) via the 
Laplace image 

= — / e''^^ ^ = ^ /e-^'/(^*) sinrdr = e-*erfi(v^) . 

27rj J Js{s + 1) Wnt J 

Br ^ ^ ' ^0 
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Fig. 1. Fractional one-half-order integral of the exponential function, 
= jV2exp(t). 

Note that erfi(x) = Jq exp^y"^) dy is the imaginary error function [6J. It is 
a real-valued, entire function defined by erfi(x) = erf(jx)/j. In contrast to 
the ordinary error function erf(x) = /g^exp(— y^) dy, the imaginary error 
function is not bounded, but the function ii{t) is bounded. It tends to zero in 
the limit t — )■ oo. Moreover, the latter is very close to the Dawson's integral 
D{x) = Jq exp^y"^) dy [7]. In this connection it should be mentioned the 
asymptotic behavior of the imaginary error function 

erfifx) = ^e^^ ( + -x~^ + -x^^ + — x^"^ + . . . 
^/^T \ 2 4 8 J 

It may be advantageous for the estimation of ii{t) with t ^ 1. To sum up 
above is the numerical simulation of ii{t) represented in Fig. [T] 



3 Decomposition 

The next important step consists in the decomposition ia{t) into two contri- 
butions. For this purpose we bow the Bromwich path of integration into the 
equivalent Hankel path. Then the loop will start from minus infinity along the 
lower side of negative real axis, encircle |s| = 1 counter-clockwise and end at 
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Fig. 2. Fractional oscillations for a = 1.8 . 

minus infinity along the upper side of tlie negative real axis. The first contri- 
bution arises from two borders of the cut negative real axis. Taking s = re^'^ 
along the upper border and s = re~^'^ along the lower border, we get 

oo 

Kit) = J e-'''V^{r)dr 





with 



Va{r) = - 



TC 



r2a 



+ 2r° cosfvra) + 1 



The part ha{t) is not completely monotonic, as the function V^(r) takes on 
both positive and negative values (see Theorem 7 in [8]). It becomes vanish- 
ingly small with t tending to infinity. The second contribution is calculated 
by means of residues. The poles Sq = exp{jn/a) and Si = exp{—jn/a) give 



J t cos{n/a) 



sm 



a 



This part clearly demonstrates an oscillatory character. Due to it, the oscilla- 
tory behavior is carried over to the function ia{t) itself. 
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It is useful to remark that the functions ia(t) resembles Cait). So for 1 < a < 2 
the function Cait) also consists of two parts, Cait) = fa(t) + Qaif)- From [Ij it 
follows that the monotonic part is 

oo 


where 

K (r) = - ^^^(^^) 

TT r2" + 2r° cos(7ra) + 1 ■ 

According to [8j, the function fa{t) is complete monotonic. The second con- 
tribution of Ca {t) takes the oscillatory form 

gjt) = -e*™^("/") COS 

a 

To have a clear idea of ia{t) and ea{t), their pictures are represented in Fig. O 



4 Equations of fractional oscillations 

The first approach to fractional oscillations had enough a formal character. 
It proceeded from the simple change of the second derivative in the harmonic 
oscillator equation to the derivative of a fractional one. In other words, such 
a equation was postulated. After solving it the fractional oscillation features 
were established. Nevertheless, there is an alternative way to fractional calcu- 
lus [9]. The probabilistic point of view shows that the derivative of fractional 
order is connected with a-stable probability distributions |T0]. Starting from 
the analysis of stochastic random processes, this permits ones to derive the 
fractional differential equations responsible for fractional oscillations. Here we 
intend to give new examples demonstrating this approach. 

Let the sign be for the juxtaposition of a function depending on t with its 
Laplace image depending on s, namely 

u{t) -T- u{s) 




oo 

J e-''u{t) dt. 
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The reader is reminded that 



+ 1 



If the Laplace image is written as 
1 1 



+ 1 s (1 + l/s°) 



the expression 

ea s + — e„ s = - 
s 



is treated by a (fractional) integral equation after the Laplace transform. Since 
6^(0) = 1, from the Laplace inversion we find 

e«(t) = 1 - J"e,(t) . (2) 



According to [TT], the definition of fractional derivative is 



taking the positive integer number m with m — 1 < a < m. In fact, D"^ is the 
m-derivative under the integral J™~". Then the function ea{t) satisfies also 
the equivalent equation in the differential form 

^"e„(t) + e„(t) = 0. (3) 



It should be pointed out, though this equation was known earlier, it was 
not derived but postulated. Using the approach, now we intend to derive an 
equation describing ia{t). 

In order to obtain the equation for ia(t), we need to consider its Laplace image 

The initial condition is ^0,(0) = 0. Thus, the corresponding equation becomes 
D'^/Xit) + J^^^Xit) = 1 . (4) 
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It is interesting to observe that following the argument of |3], the generalized 
momentum of the fractional oscillator takes the form 



where Qa is the displacement, u the circular frequency, m the generalized mass. 
To put it in another way, the momentum is expressed in terms of —ia{t)- 
Really, expanding ea{t) and ia{t) in an infinite power series about and 
differentiating the functions by with respect to t , we obtain 

whereas 

J"/'e„(t) = ia{t) , J"/'z„(t) = 1 - e„(t) . 
Consequently, Eq.(jl]) is reduced to 

D'^iait) + ia{t) = . 



The form of this equation coincides with Eq.([3]). Thus, the functions ea(t) and 
iait) represent two independent solutions of the same equation. 

Moreover, this permits one to derive the equation of fractional oscillator by 
means of the Hamilton formalism with some alterations. So the generalized 
Hamiltonian is written as 

'Ha = {pl + u'ql)/2. (5) 

It corresponds to the total energy of the fractional oscillator [3] . In the Hamil- 
ton form the motion equation of fractional oscillator becomes 



^ ^"=^ = = 

The modification of the conventional representation of the Hamilton equations 
is conditioned on a random interaction of the system of harmonic oscillators 
with environment. Each harmonic oscillator is governed by its own internal 
clock. Although its dynamics is described by the ordinary Hamiltonian equa- 
tions, p{t) and g(r) depend on the operational time. The passage from the 
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operational time to the physical time via the averaging procedure accounts for 
the interaction of the system of harmonic oscillators with environment. This 
results in this system on the whole behaving as a fractional oscillator. The 
momentum of the fractional oscillator and its displacement take the integral 
form 

oo oo 

Pait) = J p^{t, r) p(r) dr , q^{t) = J p^{t, r) g(r) dr . 





Therefore, the fractional oscillator is an ensemble average of harmonic oscil- 
lations (see details in Although the Hamiltonian ^ is not an explicit 
function of time, for non-integer values a the dynamic system is nonconserva- 
tive because of the fractional derivative of momentum. 



As is well known [12] , the Laplace images are also useful for finding an asymp- 
totic behavior of transformed functions. With this in mind we expand the 
term l/(s° + 1) in a series with respect to s. This technique works finely for 
the solutions of linear fractional differential equations [T1I13H14] . Taking the 
leading terms, we obtain the asymptotic representations as t — )■ oo, namely 

ea{t) — -, ia{t) ^ — — . 

r(l — a) r(l — a/z) 



For 1 < a < 2 the denominator r(l — a) is less than zero, whereas r(l — 
a/2) > 0. It is easy to show that the asymptotic algebraic decay of eaif) and 
iaif) is completely determined by the parts faif) and ha(t), respectively (see 
Section [3]). In the classical case a = 2 the contributions fa{t) and ha(t) are 
equal to zero exactly, because there is only a pole, and the cut on the negative 
real axis is no longer present. 



5 Zeros of fractional oscillations 



The fractional oscillations occupy a special place in the oscillatory theory and 
what is why. On the one hand, they have a decay like a relaxation. On the 
other hand they show themselves as oscillations. Asymptotically the fractional 
oscillations become vanishingly small with t oo. However, the feature has 
a deep cornerstone. Their representation in the form of two contributions 
is especially important. The decomposition manifests a competition between 
two different dependencies. The harmonic oscillation exponentially decays on 
the background of a slow algebraic relaxation. Although at the beginning the 
oscillations prevail, the algebraic decay will survive them. As a result, the 
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fractional oscillations have a finite number of damped oscillations as well as 
zeros. 



Consider the problem in greater detail. Probably, Wiman was the first who 
has described the position of the Mittag-Leffier function zeros in the complex 
plane [15] . A finite number of the zeros has been established. Later the analysis 
was specified in [T]. The function ea{t) demonstrates an odd number of zeros. 
The smallest zero lies in the interval < t < 7r/[sin(7r/a)]. The number of 
zeros strictly depends on the index a. Increasing a more and more, we get 
more the number of zeros in the function ea{t). For a = 1 the only zero is 
located in infinity. When a tends from 1 to 2, the number of zeros increases 
so that for a = 2 their number becomes equal to infinity. 

We provide a similar analysis to the function ia{t)- It is not evident that 
the function has a finite number of zeros for 1 < a < 2 . At once it should 
be noticed that this function always is zero for t = under any value of 
1 < a < 2. When t is enough large, the zeros of ia(t) are expected to be 
found approximately from the equation 



a r(l-a/2) 



neglecting the oscillation factor in the contribution qa{t)- Putting now a 
1 + e, the first-order approximation gives 

cos(7r/a) = cos[7r/(l + e)] ~ cos[7r(l — e)] ~ —1 . 



Thus, the asymptotic position T of the largest zero is defined by the relation 
2r(l/2 + 5/2) ' 



which shows that T tends to infinity as e — )■ 0. 

The function 22 (^) has infinitely many zeros. Now let the index a be 2 — 5. In 
the limit of the first-order approximation we can write 

cos(7r/a) = cos[7r/(2 - 5)] ~ cos[(7r/2)(l + 5/2)] = sin(7r5/4) ~ nS/A 



and 

r(l - a/2) = T{6/2) ~ 2/6. 
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In this case the asymptotic estimation for the largest zero T reduces to the 
equation 

g-7r5T/4 ^ ^(2 - ^) y_i+j/2 



from which one gets 

7r5T/4 - ln(2T/5) . (10) 



Since 5 — > 0, the value T tends to infinity faster than 1/5. Moreover, the terms 
7r5T/4 and ln(2T/5) are of the same order. Taking either T ~ (a/5) ln(2/(5) 
or 5 ~ 61n(T)/T, where a and h are positive constants to be determined, we 
arrive at a = b = 8/tt. The equivalent asymptotic expressions are 



Tr.— \n(lV fill 



ttS \5 

(12) 

In the limit a — >■ 2 the period of qa{t) becomes 27r. Consequently the number 
of zeros of ia{t) tends to A?" ~ T/n — > oo as 5 — > 0. 



6 Conclusions 



We have considered the fractional oscillations allied with ordinary harmonic 
ones. The main feature of the fractional oscillations is that they have a fi- 
nite number of zeros. This follows from the competition of two terms. One 
of them has an asymptotic behavior with an algebraic decay, and the other 
term contains a conventional harmonic oscillation vanishing in time because 
of an exponential decay. The second term decreases faster than this happens 
for the term with an algebraic decay. This important peculiarity of fractional 
oscillations is just reflected on the title of the paper. The analysis supports 
a key role of Mittag-Leffler functions. The fractional oscillations satisfy a lin- 
ear fractional differential equation. We have derived it, using the method of 
the Laplace transform. Our analytical solutions are confirmed by a numerical 
treatment. 
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